Intermittent polaron dynamics: Born-Oppenheimer out of equilibrium 
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We consider the non-equilibrium dynamics of a molecular level interacting with local phonon 
modes in the case of a strong polaronic shift which prevents a perturbative treatment of the problem. 
Instead, we find that in an adiabatic regime when the electronic states react faster than the phonon 
modes it is possible to provide a fully non-perturbative treatment of the phonon dynamics including 
random noise and dissipation. The result shows intermittent switching between bistable states of 
the oscillator with an effective random telegraph noise. 



The adiabatic, or Born-Oppcnhcimcr approximation is 
a key approximation in dealing with a variety of strong- 
coupling equilibrium problems in solid state, quantum 
chemistry and numerous other fields. It is based on di- 
viding an interacting system into two coupled subsys- 
tems, one of which has a much lower characteristic fre- 
quency than the other, such as splitting the description of 
a molecule into nuclear and electronic degrees of freedom. 
As a result one argues that the inelastic processes are 
suppressed by the difference in energy scales and there- 
fore the dynamics of the slower system is Hamiltonian, 
with an effective potential given by the faster system. 

In non-equilibrium situations, such as molecular tran- 
sistors or molecular junctions 0, S Bj a straightfor- 
ward application of the Born-Oppenheimer approxima- 
tion will, however, lead to incorrect or, at least, incom- 
plete physical picture. While in many cases it is rea- 
sonable to assume that the mechanical motion of the 
molecule between the leads is slow compared to the dy- 
namics of electrons hopping on and off the molecular or- 
bitals, the energy of the electronic excitations (induced 
by the applied bias voltage) may be comparable with 
the energy scales of molecule's mechanical motion, As a 
result, inelastic processes may occur, preventing a Hamil- 
tonian description of molecule's mechanical motion; in- 
stead, the effective molecular dynamics will include noise 
and damping. 

A simple model, which incorporates the essential 
physics of the problem at hand can be described by the 
following Hamiltonian: 

H M = (eo + Xx)dJd + £- + — ^— . 1 
2m 2 

Eq. represents a mechanical oscillator coupled to 
the extra charge Sd on a localized molecular orbital 
eo- We consider a case of sufficiently big A, so that 
the phonon deformation energy (or the polaronic shift) 
e p = A 2 / (2TOW 2 ,) > T, where T is the resonant level width 
defined below. The large polaronic shift will prevent us 
from working perturbatively in A. Instead, we will find 
below that it is possible to work fully non-perturbatively 
in the classical dynamics of the oscillator in an adia- 
batic limit where fluctuations about the classic trajec- 
tory can be treated perturbatively. This calculation may 
be viewed as a generalization of the Born-Oppenheimer 
approximation to this case far from equilibrium. 



The coupling of the molecule to the leads is described 
by the Hamiltonian 

H T = ^ e kac\ a c ka A a (d)c ka + c[ a dj , (2) 

k,a k,a 

where operators c ka and d describe spinless fermions in 
the left and right (a = L, R) leads as in Fig. 1 and 
and at the molecular orbital respectively. For simplic- 
ity the hybridization matrix elements A Q between the 
leads and the molecule are assumed to be independent 
of the single particle states in the leads (labelled by in- 
dex k in Eq. (J2J). It is convenient to define partial ef- 
fective hybridization widths of the molecular orbital by 
T a = 7ri/ Q A 2 , where v a are densities of states in the leads 
(the full hybridization width of the molecular orbital is 
r = Tt, + Tfl). The chemical potentials in the leads (hl 
and are assumed to be biased by an external voltage 
eV = 1-J.l — f-R- The leads are assumed to be in zero 
temperature thermostats. 

The full Hamiltonian, H — Hm + Ht, together with 
the boundary conditions fixed by the chemical potentials 
of the leads, provides a full description of the system. 
Here we are interested in the dynamics of the molecule 
as well as in the statistical characteristics of the current 
through the molecule under the strongly non-equilibrium 
conditions specified below. In particular we find that un- 
der the appropriate conditions the current through the 
system exhibits random telegraph noise. We also explic- 
itly evaluate parameters of the noise in terms of the pa- 
rameters of the problem stated above. 

Before going into the detailed analysis of the prob- 
lem we note that in the regime of strong electron-phonon 
coupling, e.g., e p 3> T, the molecule is bistable, hav- 
ing metastable states with with either and 1 charge 
occupancy of the molecular orbital. Indeed, for a neg- 
ligibly small r, the energy states of the molecule with 
an empty electronic orbital are = fuvo(n + 1/2), 
n = 0, 1, while for the occupied orbital the energies 
are E X M = e — e P + huj (n + 1/2), e.g., Eq. QJ. These 
two states are separated by an effective potential barrier 
of height ~ e p (for e p ~ eo). 

There are two possible mechanisms for the transition of 
the molecule from one state to another, either the under- 
barrier tunnelling of the molecule and or the absorption 
of energy from the electronic subsystem. Transport of 
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a single electron through the system can be accompa- 
nied by the energy transfer to the vibrational degree of 
freedom of the molecule, i.e., the oscillator, by a value 
no greater than V. Therefore, for sufficiently low volt- 
age V < wo the system is in equilibrium and the only 
available mechanism for its transition is through quan- 
tum tunnelling. 

A problem of this sort has been considered in the lit- 
erature [f|. It has been shown that the partition func- 
tion of the system can be mapped onto that of inter- 
acting Coulomb gas in 1 dimension. The properties of 
the latter are well established based on the Renormal- 
ization Group (RG) arguments [6|- The transition rate 
between the two metastable states of the molecule is pro- 
portional to the fugacity in the Coulomb gas model. For 
a symmetric case, eo = e p , the rate is proportional to 
exp(— e p /uj), which is essentially the quasiclassical tun- 
nelling exponent. An additional logarithmic renormal- 
ization of the exponent comes from the orthogonality 
of electronic states corresponding to different position- 
ing of the oscillator p| . It can be obtained from the RG 
equations [6|. The case of finite, but small bias voltage 
(V < wo), has recently been considered in Ref. 0]. 

In this paper we are primarily interested in the oppo- 
site case, V > wo ■ This corresponds to a strongly out-of- 
cquilibrium situation as the molecule can now efficiently 
absorb energy from electrons. If V is small compared to 
e p , i.e., the height of the potential barrier separating the 
two states, the system can depart from a given state only 
as result of a multi-quantum adsorption. A straightfor- 
ward perturbation expansion in powers of A clearly can 
not be applicable in this limit, but the appropriate adia- 
batic expansion will be found below. 

We wish to compute a propagator for the density ma- 
trix of the oscillator defined as 

Z osc = Tr cl [p c iT K S{-co,oo)S(co, -co)} /Tr [p c i\ . (3) 

Here <S(oo, -co) and <S(— oo, oo) are scattering operators 
for the full system, <S(oo, — oo) — exp[— i J_ oo Hdt} 1 7k 
denotes the standard time ordering along the Keldysh 
contour 8] . The density matrix of the unperturbed elec- 
trons is the direct product of the uncoupled density ma- 
trices of electron reservoirs in the leads (pl and pn) with 
an empty electron state in the resonant level (pd = deft), 
Pel = Pl ® Pr® Pd- The trace in Eq. © is taken over 
the electronic degrees of freedom. 

The oscillator effective action can be expressed in 
terms of a functional integral 
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T>x c Vx q e iSo+s * 



(4) 



where So is a bare action of the oscillator, So = 
J dt[mx c + mujQX c ]x q , and S* results from the trace over 
the electrons. One important simplification that will en- 
able us to work nonperturbatively in x c is that S* van- 
ishes identically for x q = 0, so that we can perturbatively 
expand S* in powers of x q for a given time-dependent 



x c {t). Using the cumulant expansion, the perturbative 
expansion of S* to a given order in x q can be expressed 
as a sum of connected diagrams, with Keldysh Green's 
functions which are 2x2 matrices in the space of forward 
and return indices. The vertices are x q (t) x a z , where 
a z = diag(l, — 1) in Keldysh indices. The Green's func- 
tions solve the Dyson equation G _1 = Gq 1 + Xx c (t) x 1. 

The solution of this Dyson equation is a difficult prob- 
lem. To make further progress, we exploit the smallness 
of the adiabatic parameter fl/V, where fi is a charac- 
teristic frequency of the oscillator. For V ^ the un- 
perturbed Green's function for the d electron is oscil- 
latory in time with period ~ V . Therefore the dif- 
ference in the time arguments of each x q (t) in a given 
term in S* should not significantly exceed V~ x \ define 
to to be the largest of these time arguments. Let us set 
x c (t) — x c (to) + x c (to){t — to) + .... Then, to compute a 
given term in S* , we first compute the Green's functions 
in the presence of a fixed x c (t) = x c (to) and then pertur- 
batively expand the Green's functions in derivatives of 
x c (to). One finds that a term with n derivatives of x c is 
the suppressed by (Sl/V) n . 

The x q component of x(t) essentially defines deviations 
of the oscillator's density matrix from the diagonal ma- 
trix. It can be easily checked that in an adiabatic ap- 
proximation contribution of a bubble containing n x q - 
vertices scales as (Xx q /T) n V/fl for n > 1. Therefore if 
we wish to integrate over x q in Eq. Q, we have to per- 
form an integral of a sort J dx q exp[/3x q + (V/Q,)rj(x q )]. 
In the spirit of the saddle point approximation, the fac- 
tor (V/Cl) in front again suppresses higher powers of x q 
by powers of (fl/V). 

Thus, the same parameter Cl/V controls the higher- 
derivatives of well as the higher powers of x q . 
Later, the validity of this expansion must be verified 
by a self-consistency check. We will find that in the 
region of interest, i.e., in case of two wells separated 
by a high barrier (see Fig. 1) Q ~ ojq and therefore 
the expansion is well justified far from equilibrium with 
il << V. Then, to leading order in £l/V the effec- 
tive action in Eq. involves three diagrams: the di- 
agram with a single a^-vertex, the diagram with one x q - 
vertex and one x c (to)(t — to) _v ertex, and the diagram 
with two :r g -vertices. One has to evaluate the contribu- 
tion of each of these three diagrams with the renormal- 
ized adiabatic Green's function of the d electron, defined 
by G~ x — Gq 1 + \x c x 1. The off-diagonal components 

of G, i.e., -i(df(t)dl(t )) and -i(d r (t)d j f (t )) can be ex- 
plicitly expressed in terms of their Fourier transforms, 
(l/2m) J du exp[-iu(t - t )]Gij[u, x c (t )], i,j = f,r: 



Gf r [u>, x c ] = 2i 



T L 0{p L 



F R 0(pR - LO) 



(w-e -A.T c ) 2 + r 2 



(5a) 



n r i . r L e(w - p L ) + T R e(uj - mi) 

G ^' Xc] = (co-e -\x c r+T2 ' (5b) 

In terms of the adiabatic Green's function G the semi- 
classical effective action for the oscillator can be written 
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cLS S-_ 



dt[iF(x c )x q + iA(x c )x c x q + (\/2)D(x c )x 2 



where 



F{x c ) = [ dwGf r [u, x c ] , 
2m J 

X 2 f d 
A(x c ) = — / du>Gf r [uj,x c } — G r f[w, x c 

A 2 f 

D(x c ) = — I dujGf r [uj, X c ]G r f[ui, x c ] . 
2n J 



(6) 

(7a) 
(7b) 
(7c) 



The first term in Eq. @ is a classical force exerted 
by electrons onto the oscillator. Therefore the effec- 
tive potential for the oscillator is U(x) — mwgi 2 /2 + 
f x dx'F(x'). In the bistable regime, e p = X 2 /(2mujQ) > 
T, < e < c p , and not too high voltage, lu < V < E p , 
U{x) develops two minima separated by a barrier. The 
minima correspond to the states of the molecule with oc- 
cupied and empty orbital. Though the approach devel- 
oped in this paper allows to deal with a pretty arbitrary 
situation as long as Q, <C V, here we will concentrate on 
the study of the interesting bistable case. For mathemat- 
ical simplicity let us assume that the hybridizations 
and Tr with leads are equal. Then the two minima are 
located at xq — eo/X and at x\ = (eo — 2e p )/X, with en- 
ergies E = —e^/Aep and E x — ~ (e — 2e p ) 2 /4e p , see Fig. 
1, while maximum is at x — 0. In the limit E p 3> F the 
wells are parabolic with frequencies of small oscillation 
are equal to loq- 
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FIG. 1: Symbolic representation of system's energetics: the 
potential curve of the mechanical degree of freedom versus 
the electronic energy levels. 

The last two terms in S* in Eq. © correspond to dis- 
sipation and fluctuation of the oscillator due to electrons. 
It is convenient to decouple the last term in Eq. JfjJ by 
introducing the Gaussian noise field £(t) in the functional 
integral Eq. |0J by means of Hubbard-Stratanovich trans- 
formation and integrate out the x q field. As a result the 



partition function 2 osc can be explicitly related to the 
distribution function for a classical particle, whose dy- 
namics is governed by a Langevin equation: 



T>£ / T>x c exp 



dt 



eit) 



2D(x c ) 



(8) 



S t [mx c - F(x c ) + A(x c )x c - £(t)] 



where the Gaussian noise £(t) is white, i.e., (£(£)£(£')) = 
D(x)8(t~t') 1 and St in Eq. JSJl is a functional (^-function. 
From now on we will drop the unnecessary c-subscript. 

The damping of the oscillator defined by an effective 
friction coefficient "f(x) = A(x)/m, where A is given 
by Eq. (j7b|l . is strongest in the vicinity of the barrier, 
7 barrier ~ £p(wo/r). Away from the barrier it decays 
as 7(2;) ~ x -4 and the system is underdamped every- 
where except for a narrow region of width T/X at the 
barrier, where it approaches a separatrix and thus its 
period diverges. Therefore its natural to consider mo- 
tion of the oscillator as a diffusion in energy space. The 
rate of change of the mean energy of the oscillator over 
a period depends of both dissipation and the fluctuation 



as {dE/dt) T 



-(A(x)x 2 ) T + {£(t)x) T . The friction- 



dependent term can be expressed in terms of the in- 
stantaneous energy of the oscillator as —1/T§ dtA(x)x 2 , 
where integral is taken over a period of motion T of the 
oscillator. Since the motion of the oscillator during one 
period is nearly deterministic (for small fluctuation and 
dissipation), the noise introduced by the (£(i)x)T term 
is also Gaussian. The first moment of this noise is given 
by 1/(2T) § dtD(x) (the Ito term), while the second mo- 
ment is 1/T § dtD(x)x 2 . The Ito term offsets the damp- 
ing term so that, due to fluctuations, the system tends 
to occupy states close to, but not exactly at, the bot- 
tom of the potential curve U. The Langevin equation for 
system's energy can be cast in a form: 



dE 
~dt 



= -7(£) + &(*). 



^E{t)iE{t')) = D{E)8{t-t'), 



(9) 



with the damping and the diffusion coefficients in the 
energy space: 



1(E) = ^Jdx 



A{x)x 



D(x) 



2,i: 



D{E) = - I dxD{x)x, 



(10a) 
(10b) 



where x = W2[E — U(x)]/m, the period of motion 
T = 2 j dx y/m/2 [E — U(x)], and the integration runs 
between the turning points defined by U(x) = E. 

Eq. © describes diffusion of the oscillator's energy in 
the nonlinear potential U. The Fokker-Plank equation 
which corresponds to dynamics in Eq. (JHJ reads: 



dP 
~dt 



d_ 
~dE 



1 d 2 



2dE 2 



[D(E)P] 



(11) 
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In case of the two- well potential there are three connected 
regions: two in each well and one above the barrier. Solu- 
tions of the Fokker-Plank equation in each region are re- 
lated through a boundary conditions Pq(E = 0)+Pi(E = 
0) = Pabove(£ = 0), and P (E = 0) = P X (E = 0), 
which state that when the system reaches the separatrix, 
it is equally likely to be in both wells. In equilibrium, 
i.e., in the absence of the net flux, the stationary distri- 
bution for the oscillator can be found from an equation 
7 P st - -(l/2)d(DP st )/dE. 

The current operator can be defined as I = (N^ ~ 
Nr)/2, where N a = Y^k c ka c ka, a — L,R and N a = 
i[N a ,H]. In the adiabatic limit the current expectation 
value can be expressed in terms of the stationary dis- 
tribution of the oscillator as J dEP st (E)(I(E)). Here 
(1(E)) = l/T § dt(I(x)), where (I(x)) is the expectation 
value of the current evaluated for a fixed position x of 
the oscillator (brackets stand for averaging over the elec- 
tronic degrees of freedom only) . Evaluation of the above 
integrals is easy for e p ^> max(U, T). Then the oscillator 
tends to be near the bottom of one of the two potential 
wells with negligible probability of being near the top of 
the barrier, leading to 

(I)=hn F (e Q )+I 2 [l-n* F (e )], 
*, 1 (12) 

" F(eo) ~ i+cxp[4( eo - £p )/ir 

where h = r 2 V/[27r(e - 2e p ) 2 ] and I 2 = r 2 U/[27re 2 ] 
are the currents through the molecule for the oscillator 
positioned at the minima of each well. The coefficient 
n* F has the same form as the Fermi distribution function 
with temperature V/4 and chemical potential e p . 

In this case, we find that each well has the same effec- 
tive temperature U/4, and thus the occupation probabil- 
ity of a given well in Eq. (|12[1 is given by a Boltzmann 
weight. In general, each well may have a different ef- 
fective temperature; in that case, the oscillator may be 
more likely to be found in the higher energy state if the 
effective temperature of that well is lower. 

The current-current correlation function requires eval- 
uating the expectation value, averaged over electronic 
states, of (I(E,ti)I(E,t 2 )). In the adiabatic limit the 
connected part of this expectation value can be expressed 



in terms of the stationary distribution of the oscillator as 
JdEP st (E)(I 2 (E)) con , where I 2 (E) = l/T § dt(I 2 (x)) . 
This corresponds to the shot noise ~ 2e(I). 
The disconnected part is 

J dE 1 dE 2 P tI (E 1 ,t 1 ;E 2 ,t 2 ){I(E 1 ,t 1 )){I(E 2 ,t 2 )), (13) 

where Pf[(E\,t\; E 2 ,t 2 ) is the 2-point distribution func- 
tion for the oscillator. This part contains a contribution 
due to the oscillations in each well, and due to the tran- 
sitions between the wells. For a sufficiently high bar- 
rier (E p 3> V), the transitions (jumps) are rare events 
(note that V is playing role of temperature in the wells), 
and therefore the oscillations and the jumps occur on 
different time scales. The contribution of oscillatory mo- 
tion is pea ked at lo = loq, with a width of the order of 
7(^min)lly|- The contribution of jumps into the current- 
current correlator in Eq. fLlf is J^ij A? ( f i — t 2 )IiIjP s tj, 
where i,j = 1,2 label the states in each well, and P s a 
are the probabilities for the oscillator to be found in each 
well, P s ti = J n-dEP s t,(E). The transition probabilities 
Pij between the two wells are exponentially decaying on 
time-scales Tj, where T\ are the mean first passage times 
for each well 01 • I n this two-state approximation the 
correlation function corresponds to that of a telegraph 
noise with spectrum S^f = S tcl (0)/(l + w 2 t 2 ), where 
S tc \0) = 4rr 1 r 2 - 1 r 3 (7i - I 2 ) 2 , and r" 1 = rf 1 + r^ 1 . 
The mean first passage times Tj can be evaluated from 
the Fokker-Plank Eq. (|llfl based on standard techniques, 
i.e., by solving an equation adjoint to Eq. Ill 1H [12|. One 
finds: 

Ti ={-y) J»cxp(- — ), (14) 

where Ei = e%/4E p and E 2 = (e - 2e p ) 2 /4E p , e.g., Fig. 
1. The telegraph switching is thus highly sensitive to 
the bias voltage and can be controlled by both the bias 
voltage and the gate voltage through the variation of the 
molecular orbital energy (eo) relative to the conduction 
bands of electrons in the leads. 
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